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Abstract 

Time-dependent CP asymmetries in the decay rates of the singly Cabibbo-suppressed 
decays D° -> K~ K + and D° —» 7 r 7r + are measured in pp collision data corresponding 
to an integrated luminosity of 3.0 fb -1 collected by the LHCb experiment. The 
D° mesons are produced in semileptonic 6-hadron decays, where the charge of 
the accompanying muon is used to determine the initial state as D° or D°. The 
asymmetries in effective lifetimes between D° and D° decays, which are sensitive to 
indirect CP violation, are determined to be 

A r {K-K + ) = (-0.134 ± 0.077 , 

A r {ir~iT + ) = (-0.092 ± 0.145 , 

where the first uncertainties are statistical and the second systematic. This result is 
in agreement with previous measurements and with the hypothesis of no indirect 
CP violation in D° decays. 
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1 Introduction 


In neutral meson systems, mixing may occur between the particle and anti-particle states. 
This mixing is very small in the charm-meson ( D °) system. Experimentally, a small, 
non-zero D°-D° mixing is now firmly established by several experiments [l]^6] , where the 
average of these measurements excludes zero mixing at more than 11 standard deviations [7|. 
This opens up the possibility to search for a breaking of the charge-parity (CP) symmetry 
occurring in the D°-D° mixing alone or in the interference between the mixing and decay 
amplitudes. This is called indirect CP violation and the corresponding asymmetry is 
predicted to be 0( 10 -4 ) (8,91, but can be enhanced in theories beyond the Standard 
Model 10 . Indirect CP violation can be measured in decays to CP eigenstates such as 


the singly Cabibbo-suppressed decays D° —> K~ K + and D° —> (the inclusion of 

charge-conjugate processes is implied hereafter) from the asymmetry between the effective 
D° and D° lifetimes, Ay . The effective lifetime is the lifetime obtained from a single 
exponential fit to the decay-time distribution. Several measurements of Ay exist (T il 12 


The most precise determination was made by LHCb with data corresponding to 1.0 fb~ 
of integrated luminosity, resulting in A r (K~ K + ) = (—0.035 ± 0.062 ± 0.012)%, and 
%r(7r _ 7r + ) = (0.033 ± 0.106 ± 0.014)% |11|. When indirect CP violation is assumed to 
be the same in the two modes, the world average becomes Ay = (—0.014 ± 0.052)% [7|. 
In all previous measurements of Ay, the initial flavour of the neutral charm meson (i.e., 
whether it was a D° or D° state) was determined (tagged) by the charge of the pion in 
a D* + —>■ D° 7r + decay. In this paper, the time-dependent CP asymmetry is measured 
in D° decays originating from semileptonic 6-hadron decays, where the charge of the 
accompanying muon is used to tag the flavour of the D° meson. These samples are 
dominated by B~ —)■ and B° —> D°/i~V^X decays, where X denotes other 

particle(s) possibly produced in the decay. The same data samples as for the measurement 
of time-integrated CP asymmetries 13 are used. 


2 Formalism and method 

The time-dependent CP asymmetry for a neutral D meson decaying to a CP eigenstate, 
/, is defined as 


_ r(D 0 ^f;t)-r(D°^f-,t) 
^(D°^f■t) + ^(Dv^f■t) , 


(i) 


where T(D° /; t) and r(Z7° —> /; t ) are the time-dependent partial widths of initial D° 
and D° mesons to final state /. The CP asymmetry can be approximated as 
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where A^p is the direct CP asymmetry and r is the D° lifetime. The linear decay-time 
dependence is determined by A?, which is formally defined as 


At 


r D o - r 50 

f d° + r pa 


(3) 


where T is the effective partial decay rate of an initial D° or D° state to the CP eigenstate. 
Furthermore, At can be approximated in terms of the D°-D° mixing parameters, x and y, 
as (l5) 

A r ~ (%™% /2 — A^p) y cos 0 — x sin 0 , (4) 


where = \q/p\ 2 — 1 describes CP violation in D°-D a mixing, with q and p the 


coefficients of the transformation from the flavour basis to the mass basis, |.D li2 ) = 
p\D°) ± q\D °). The weak phase 0 describes CP violation in the interference between 
mixing and decay, and is specific to the decay mode. Finally, At receives a contribution 


from direct CP violation as well 16 


The raw asymmetry is affected by the different detection efficiencies for positive and 
negative muons, and the different production rates of D° and D° mesons. These effects 
introduce a shift to the constant term in Eq. (|2j, but have a negligible effect on the 
measurement of At (see Sect. [6J. The decay D° —> K~n + , also flavour-tagged by the 
muon from a semileptonic 6-hadron decay, is used as a control channel. Since this is a 
Cabibbo-favoured decay mode, direct CP violation is expected to be negligible. More 
importantly, any indirect CP violation is heavily suppressed as the contribution from 
doubly Cabibbo-suppressed D° —» K + tt~ decays is small. 


3 Detector and simulation 

The LHCb detector (l7p~8;j is a single-arm forward spectrometer covering the pseudorapidity 
range 2 < 77 < 5, designed for the study of particles containing b or c quarks. The detector 
includes a high-precision tracking system consisting of a silicon-strip vertex detector 
surrounding the pp interaction region, a large-area silicon-strip detector located upstream 
of a dipole magnet with a bending power of about 4Tnr, and three stations of silicon-strip 
detectors and straw drift tubes placed downstream of the magnet. The polarity of the 
magnetic held is regularly reversed during data taking. The tracking system provides a 
measurement of momentum, p, with a relative uncertainty that varies from 0.5% at low 
momentum to 1.0% at 200GeV/c. The minimum distance of a track to a primary vertex, 
the impact parameter, is measured with a resolution of (15 + 29/pp) pm, where p t is 
the component of the momentum transverse to the beam, in GeV/c. Different types of 
charged hadrons are distinguished using information from two ring-imaging Cherenkov 
detectors. Photon, electron and hadron candidates are identified by a calorimeter system 
consisting of scintillating-pad and preshower detectors, an electromagnetic calorimeter and 
a hadronic calorimeter. Muons are identified by a system composed of alternating layers 
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of iron and multiwire proportional chambers, situated behind the hadronic calorimeter. 
The trigger 19 consists of a hardware stage, based on information from the calorimeter 


and muon systems, followed by a software stage, which applies a full event reconstruction. 

In the simulation, pp collisions are generated using Pythia 
configuration (2l]. Decays of hadronic particles are described by EvtGen 
final-state radiation is generated using PHOTOS 23 . The interaction of the generated 


with a specific LHCb 
in which 
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particles with the detector, and its response, are implemented using the Geant4 toolkit 
as described in Ref. 
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4 Data set and selection 


This analysis uses a data set corresponding to an integrated luminosity of 3.0 fb 1 . The 


data were taken at two different pp centre-of-mass energies: 7TeV in 2011 (1.0 fb x ) and 


8TeV in 2012 (2.0fb^ ). The data sets recorded with each dipole magnet polarity are 
roughly equal in size. 

At the hardware trigger stage, the events are triggered by the presence of the muon 
candidate in the muon system. This requires the muon p t to be greater than 1.64GeV/c 
(1.76GeV/c) for the 2011 (2012) data. At the software trigger stage, one of the final-state 
particles is required to have enough momentum and be significantly displaced from any 
primary pp vertex. In addition, the candidates must be selected by a single-muon trigger 
or by a topological trigger that requires the muon and one or two of the D° daughters to 


be consistent with the topology of 5-hadron decays 19 . 

To further suppress background, the D° daughters are required to have pt > 300 MeV/c. 
All final-state particles are required to have a large impact parameter and be well identified 
by the particle identification systems. The impact parameter requirement on the muon 
reduces the contribution from D° mesons produced directly in the pp interaction to below 
2%. The scalar p t sum of the D° daughters should be larger than 1.4GeV/c, and the 
Px of the D° candidate should be larger than 0.5GeV/c. The two tracks from the D° 
candidate and the D°p combination are required to form good vertices and the latter 
vertex should be closer to the primary vertex than the D° vertex. The D° decay time 
is determined from the distance between these two vertices, and the reconstructed D° 
momentum. The invariant mass of the D°p combination is required to be between 2.5 and 
5.0GeV/c 2 , where the upper bound suppresses hadronic 5-hadron decays into three-body 
final states. Backgrounds from inclusive 5-hadron decays into charmonium are suppressed 
by vetoing candidates where the invariant mass of the muon and the oppositely charged D° 
daughter, misidentihed as a muon, is consistent with the J/0 or ip{2S) mass. Additionally, 
the invariant mass of the muon and same-charge D° daughter, under the muon mass 
hypothesis, is required to be larger than 240 MeV/c 2 to remove events where a single 
charged particle is reconstructed as two separate tracks. For most selection requirements, 
the efficiency is roughly independent of the D° decay time, giving efficiency variations of 
0(1%). The largest dependence on the decay time comes from the topological trigger, 
which introduces an efficiency profile that decreases with D° decay time, resulting in about 
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20% relative efficiency loss at large decay times. 


5 Determination of Ay 


The mass distributions for the selected D° —> K~K + , D° —> 7r - 7r + and D° —> K~n + 
candidates are shown in Fig. [lj The numbers of signal candidates are determined from 
unbinned extended maximum-likelihood fits in the range 1810 to 1920MeV/c 2 . The signal 
for all three decay modes is modelled by a sum of three Gaussian functions. The first two 
have the same mean, but independent widths; the third is used to describe a small radiative 
tail, and has a lower mean and larger width. The effective width of the signal ranges 
from 7.1MeV/c 2 for D° —> K~K + candidates to 9.3MeV/c 2 for D° —> 7r~n + candidates. 
As the final states Ii~K + and 7r - 7r + are charge symmetric, the shape parameters for the 
signal are the same for both D° and D° candidates. The combinatorial background is 
modelled by an exponential function. In the 7r“7r + invariant mass distribution, a reflection 
from D° —> K~tt + decays is visible in the region below 1820MeV/c 2 . This background 
component is modelled by a single Gaussian function and the fit range is extended from 
1795 to 1930MeV/c 2 . The shape parameters and overall normalisation of the background 
components are allowed to differ between D° and D° candidates. The numbers of signal 
candidates obtained from these global fits are 2.34 x 10 6 for D° —> Ii~K + , 0.79 x 10 6 for 
D° —>• 7r - 7r + and 11.31 x 10 6 for D° —> K n + decays. The latter number corresponds to 
only half of the available D° —>• K~n + candidates, randomly selected, to reduce the sample 


size. 


The raw CP asymmetry is determined from fits to the mass distributions in 50 bins 
of the D° decay time. The fits are performed simultaneously for D° and D° candidates 
and the asymmetry is determined for each decay-time bin. The shape parameters and 
relative normalisation for the third Gaussian function and for the D° —> K~ n + reflection 
background are fixed from the global fit. All other parameters are allowed to vary in these 
fits. In particular, since both the amount and the composition of background depend on the 
decay time, the background parameters are free to vary in each decay-time bin. For decay 
times larger than 1 ps the relative contribution from combinatorial background increases. 
This is due to the exponential decrease of the signal and a less steep dependence of the 
combinatorial background on the decay time. The mass distribution in each decay-time 
bin is well described by the model. 

Events at large D° decay times have a larger sensitivity to Ap compared to events 
at small decay times, which is balanced by the fewer signal candidates at large decay 
times. The binning in D° decay time is chosen such that every bin gives roughly the 
same statistical contribution to Ap. The value of Ap is determined from a y 2 fit to the 
time-dependent asymmetry of Eq. ([2]). The value of Ap and the offset in the asymmetry 


are allowed to vary in the fit, while the D° lifetime is fixed to r = 410.1 fs 26 . Due to the 


exponential decay-time distribution, the average time in each bin is not in the centre of the 
bin. Therefore, the background-subtracted |27| average decay time is determined in each 
bin and used in the fit. This fit procedure gives unbiased results and correct uncertainties, 
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Figure 1: Invariant mass distributions for (a) D° —> K~K + , (b) D° —> tt~tt + and (c) D° —> K~n + 
candidates. The results of the fits are overlaid. Underneath each plot the pull in each mass bin 
is shown, where the pull is defined as the difference between the data point and total fit, divided 
by the corresponding uncertainty. 


as is verified by simulating many experiments with large samples. 

The measured asymmetries in bins of decay time are shown in Fig. [2j including the 
result of the time-dependent fit. The results in the three decay channels are 

A r (K~K+) = (-0.134 ±0.077)% , 

A r (vr _ 7r + ) = (-0.092 ± 0.145)% , 

Ar(K~TT + ) = ( 0.009 ± 0.032)% , 

where the uncertainties are statistical only. The values for Ap are compatible with the 
assumption of no indirect CP violation. The fits have good p-values of 54.3% (D° —y 
K~K + ), 30.8% (D° —y and 14.5% (D° —y K~ir + ). The measured values for the raw 

time-integrated asymmetries, which are sensitive to direct CP violation, agree with those 
reported in Ref. [l3j. 
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Figure 2: Raw CP asymmetry as function of D° decay time for (a) D° —» K~K + , (b) D° —> ir~ir + 
and (c) D° —» I\~n + candidates. The results of the x 2 are shown as blue, solid lines with the 
±1 standard-deviation (<r) bands indicated by the dashed lines. The green, dashed lines indicate 
one D° lifetime (r = 410.1 fs). Underneath each plot the pull in each time bin is shown. 
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Table 1: Contributions to the systematic uncertainty of Ay(K K + ) and Ar(vr 7r + ). The constant 
and multiplicative scale uncertainties are given separately. 


Source of uncertainty 

D° in¬ 
constant 

-K+ 

scale 

D°—)• t r 
constant 

7T + 

scale 

Mistag probability 

0.006% 

0.05 

0.008% 

0.05 

Mistag asymmetry 

0.016% 


0.016% 


Time-dependent efficiency 

0 .010% 


0 .010% 


Detection and production asymmetries 

0 .010% 


0 .010% 


D° mass fit model 

0 .011% 


0.007% 


D° decay-time resolution 


0.09 


0.07 

B°-B° mixing 

0.007% 


0.007% 


Quadratic sum 

0.026% 

0.10 

0.025% 

0.09 


6 Systematic uncertainties and consistency checks 

The contributions to the systematic uncertainty on Ay are listed in Table [l] The largest 
contribution is due to the background coining from random combinations of muons and 
D° mesons. When the muon has the wrong charge compared to the real D° flavour, this 
is called a mistag. The mistag probability (cu) dilutes the observed asymmetry by a factor 
(1 — 2cn). This mistag probability is measured using D° —> K~n + decays, exploiting the 
fact that the final state determines the flavour of the D° meson, except for an expected 
time-dependent wrong-sign fraction due to D°-D° mixing and doubly Cabibbo-suppressed 
decays. The mistag probability before correcting for wrong-sign decays is shown in Fig. [3j 
After subtracting the (time-dependent) wrong-sign ratio |3|, the mistag probability as 
function of D° decay time is obtained. The mistag probability is small, with an average 
around 1%, but it is steeply increasing, reaching 5% at five D° lifetimes. This is due to the 
increase of the background fraction from real D° mesons from 6-hadron decays combined 
with a muon from the opposite-side 6-hadron decay. This random-muon background is 
reconstructed with an apparently longer lifetime. The time-dependent mistag probability 
is parameterised by an exponential function, which is used to determine the shift in Ay . 
The systematic uncertainty from this time-dependent mistag probability is 0.006% for 
the D° —>- K~K + and 0.008% for the D° —7r - 7r + decay mode, with a supplementary, 
multiplicative scale uncertainty of 0.05 for both decay modes. 

The mistag probabilities can potentially differ between positive and negative muons. 
Such a mistag asymmetry would give a direct contribution to the observed asymmetry. 
The slope of the mistag asymmetry is also obtained from D 0 —> K~7 r + decays. This slope 
is consistent with no time dependence, and its statistical uncertainty (0.016%) is included 
in the systematic uncertainty on Ay . 

The selection of signal candidates, in particular the topological software trigger, is 
known to introduce a bias in the observed lifetime. Such a bias could be charge dependent, 
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Figure 3: Mistag probability, before subtracting the contribution from wrong-sign (WS) decays, 
determined with D° —> K~ir + candidates. The result of the fit to the data points with an 
exponential function is overlaid (solid, blue line). The red, dashed line indicates the expected 
mistag contribution from WS decays. 


thus biasing the measurement of A p. It is studied with the D° — > K~tt + sample and a 
sample of D~ —> K + n~7i~ decays from semileptonic 6-hadron decays. No asymmetry of 
the topological triggers in single-muon-triggered events is found within an uncertainty of 
0.010%. This number is propagated as a systematic uncertainty. 

The detection and production asymmetries introduce a constant offset in the raw 
time-dependent asymmetries. Since these asymmetries depend on the muon or 6-hadron 
momentum, they can also introduce a time dependence in case the momentum spectrum 
varies between decay-time bins. This effect is tested by fitting the time-dependent 
asymmetry after weighting the events so that all decay-time bins have the same D° 
or muon momentum distribution. The observed shifts in Ay are within the statistical 
variations. The shift (0.010%) observed in the larger D° — )■ K~tt + sample, which has the 
same production asymmetry and larger detection asymmetry, is taken as a measure of the 
systematic uncertainty. 

An inaccurate model of the mass distribution can introduce a bias in Ap. The effect 
on the observed asymmetries is studied by applying different models in the fits to the 
invariant mass distributions. For the signal, a sum of two Gaussian functions with and 
without an exponential tail, and for the background a first and a second-order polynomial 
are tested. The maximum variation from the default fit for each decay mode (0.011% for 
D° —> K~K + \ 0.007% for D° —> 7r - 7r + ) is taken as a systematic uncertainty on Ay- 

The D° decay-time resolution affects the observed time scale, and therefore changes 
the measured value of Ay . For each decay mode, the resolution function is obtained from 
the simulation, which shows that for the majority of the signal (90%) the decay time is 
measured with an RMS of about 103 fs. The remaining candidates (10%) are measured 
with an RMS of about 312 fs. The theoretical decay rates are convolved with the resolution 
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Due 


functions in a large number of simulated experiments. The effect of the time resolution 
scales linearly with the size of Ap. The corresponding scale uncertainty on Ap is 0.09 for 
the D° —» K~K + decay mode and 0.07 for the D° —)• 7r _ 7r + decay mode. Decays where the 
muon gives the correct tag but the decay time is biased, e.g., when the muon originates 
from a r lepton in the semileptonic 6-hadron decay, are studied and found to be negligible. 

About 40% of the muon-tagged D° decays originate from neutral B mesons 28 
to B°-B° mixing the observed production asymmetry depends on the B° decay time 
A correlation between the B° and D° decay times may result in a shift in the measured 
value of Ap. The effect of this correlation, determined from simulation, together with a 
1% B° production asymmetry 29.30 , is estimated to be a shift of 0.007% in the observed 
value of Ap. This is taken as systematic uncertainty. 

Possible shifts in Ap coming from the 1.5 fs uncertainty on the world-average D° 
lifetime 26 , from the uncertainty on the momentum scale and detector length scale 31,32 


29 


and from potential biases in the fit method are negligible. 

The scale uncertainty (cf. Table [l| gives a small contribution to the overall systematic 
uncertainty, which depends on the true value of Ap. In order to present a single systematic 


uncertainty, the effect of the scale uncertainty is evaluated with a Neyman construction 33 


For each true value of Ap, the absolute size of the scale uncertainty is known and added in 
quadrature to the constant uncertainty. In this way, a confidence belt of observed values 
versus true values is constructed. This procedure gives a slightly asymmetric systematic 
uncertainty, which is +o;o 34 % for the D° —> K~K + decay channel and 1 q!o 33% f° r the 
D° — > 7i~7r + decay channel. Except for the contribution from the mass fit model, all 
contributions to the systematic uncertainty are fully correlated, resulting in an overall 
correlation coefficient of 89% between the systematic uncertainties of Ap(K~K + ) and 
Ap(7T~7T + ). 

Additional checks have been performed to determine potential sensitivity of the mea¬ 
surements on the data-taking conditions, detector configuration, and analysis procedure. 
Changing to a finer decay-time binning yields compatible results. Potential effects on 
the measurement of Ap coming from detection asymmetries are expected to appear when 
dividing the data set by magnet polarity and data-taking period. Detection asymmetries 
originating from a left-right asymmetric detector change sign when reversing the magnet 
polarity. Similarly, during the two data-taking periods, detection asymmetries and produc¬ 
tion asymmetries might have changed due to different running conditions. As shown in 
Fig-S there is no significant variation of Ap across various configurations. Also splitting 
the data set according to the number of primary vertices or in bins of the B decay time 
does not show any deviation in the measured values of Ar- 


7 Conclusions 

The time-dependent CP asymmetries in D° —> K~ K + and D° —> 7i~7r + decays are measured 
using muon-tagged D° mesons originating from semileptonic 6-hadron decays in the 3.0 fb _1 
data set collected with the LHCb detector in 2011 and 2012. The asymmetries in the 


9 












Mag. up 2011 


Mag. up 2011 
Mag. down 2011 
Mag. up 2012 
Mag. down 2012 
All 2011 
All 2012 

-1 - 0.5 0 0.5 1 -1 - 0.5 0 0.5 1 

A r [%] A r [%] 

Mag. up 2011 
Mag. down 2011 
Mag. up 2012 
Mag. down 2012 
All 2011 
All 2012 

-1 - 0.5 0 0.5 1 

A r [%] 

Figure 4: Measured values of Ap for different magnet polarities and data-taking periods for 
(a) D° —> K~I\ + , (b) D° —> and (c) D° —> K~ir + decays. The vertical line and error 

band indicate the average Ap obtained from the combined data set. The error bars indicate the 
statistical uncertainty only. 
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effective lifetimes are measured to be 


A r (K~K + ) = (-0.134 ± 0.077 1°;°^)% , 
Ap(vr-7r + ) = (-0.092 ± 0.145 1^|)% , 


where the first uncertainty is statistical and the second systematic. Assuming that indirect 
CP violation in D° decays is universal 
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and accounting for the correlation in the 
systematic uncertainties, the average of the two measurements becomes Ap = (—0.125 ± 
0.073)%. The results in this paper are uncorrelated with the time-integrated asymmetries 
reported in Ref. 13 . The results are consistent with other Ap measurements (1,11,12 , and 
independent of the Ap measurements (TT| from LHCb using D° mesons from D* + —> D°7i + 
decays. They are consistent with the hypothesis of no indirect CP violation in D° —> K~K + 
and D°—t 7i~n + decays. 
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